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Abstract 

Analysis  of  overall  energy  conversion  in  laser  propulsion  is  reported.  Experimental  studies  of  a  laboratory 
scale  propulsion  device  that  absorbs  laser  energy  and  converts  that  energy  to  propellant  kinetic  energy  were  carried 
out.  The  Myrabo  Laser  Lightcraft  (MLL),  propelled  by  laser  heated  air,  was  studied.  The  MLL  incorporates  an 
inverted  parabolic  reflector  that  focuses  laser  energy  into  a  toroidal  volume  where  it  is  absorbed  by  a  unit  of 
propellant  mass  that  is  subsequently  expanded  in  the  geometry  of  the  plug  nozzle  aerospike.  Thermodynamics 
predicted  that  the  upper  limit  of  the  efficiency  of  conversion  of  the  internal  energy  of  laser  heated  air  to  jet  kinetic 
energy,  a,  is  ~  0.30  for  EQUILIBRIUM  expansion  to  1  bar  pressure.  The  analysis  captures  the  equation  of  state  of 
partially  ionized  air  under  conditions  of  chemical  equilibrium.  This  upper  limit  a  is  nearly  independent  of  the 
specific  internal  energy  between  1  and  1 00  MJ/kg,  or  temperature  from  2000  to  24000  K  at  density  of  1 . 1 8  kg/mJ. 
The  upper  limit  efficiency'  for  optimum  FROZEN  expansion  of  laser  heated  air  is  a  =  0.27.  With  heating  of  air  at 
its  Mach  5  stagnation  density  (5.9  kg/mJ  as  compared  to  STP  air  density  of  1 .1 8  kg/m')  these  efficiencies  increase  to 
about  0.55  (equilibrium)  and  0.45  (frozen).  Optimum  blowdown  from  1.18  kg/mJ  to  1  bar  occurs  with  expansion 
ratios  ~  1.5  to  4  as  internal  energy  decreases  from  1  to  100  MJ/kg.  Heating  of  Mach  5  air  at  stagnation  density 
requires  larger  expansion  ratios,  8  to  32,  for  optimum  expansion  to  1  bar.  Expansion  of  laser  ablated  Delrin 
propellant  appears  to  convert  the  absorbed  laser  energy  more  efficiently  to  jet  kinetic  energy  because  the  effective 
density  of  the  ablated  gaseous  Delrin  is  significantly  greater  than  that  of  STP  air. 

NOMENCLATURE  fin  order  of  usef 

Ef  kinetic  energy  of  vehicle  at  end  of  mission, 
nif  mass  of  vehicle  at  end  of  mission.  • 

vf  velocity  of  vehicle  at  end  of  mission  in  inertial  frame  of  reference. 

t|  efficiency  of  conversion  of  propellant  kinetic  energy  to  vehicle  kinetic  energy, 
a  efficiency  of  conversion  of  propellant  internal  energy  to  propellant  kinetic  energy. 

P  efficiency  of  absorption  of  laser  energy  by  propellant, 

y  efficiency  of  transmission  of  laser  energy  through  atmosphere  to  vehicle. 

El  laser  energy  per  laser  pulse. 

V;  initial  velocity  of  rocket  in  inertial  frame,  m/s. 

vc  exit  velocity  of  propellant  relative  to  the  rocket,  in  the  rocket  frame  of  reference,  m/s. 
m;  initial  mass  of  rocket,  kg. 

f  mass  fraction  for  a  rocket  mission,  f  =  mj/m;. 

F  Force  or  thrust  of  rocket,  N  =  kg  m/s2, 
m  mass  of  rocket,  kg. 

t  time,  s 

dm/dt  incremental  mass  change  of  rocket,  instantaneous  propellant  mass  flow  rate,  kg/s. 
p  propellant  density,  kg/m3. 

<ve>  mass  weighted  average  exit  velocity  of  blowdown  expansion  in  rocket  frame  of  reference,  m/s. 

v  velocity  of  rocket  in  inertial  frame  of  reference,  m/s. 

dv/dt  incremental  velocity  change  of  rocket  in  the  inertial  frame  of  reference,  m/s2. 

x  ratio  of  vr -Vi  to  ve. 

y  ratio  of  vs  to  ve. 

Ep  kinetic  energy  of  propellant,  J. 

<ve2>  mass  weighted  average  squared  exit  velocity  of  propellant  in  a  blowdown  expansion,  m2/s2. 
mp  mass  of  propellant,  kg. 


I  impulse,  Ns  =  kg  m/s. 

C  coupling  coefficient,  Ns/J.  . 

<I>  ratio  of  square  of  mass  weighted  average  exit  velocity  to  mass  weighted  root  mean  square  exit  velocity, 
<ve>2/<v,>. 

Uc  specific  internal  energy  of  laser  heated  propellant,  J/kg. 

u°  specific  internal  energy  of  propellant  before  laser  heating,  u0lir  at  STP  =  -  9.0  x  104  J/kg. 

u„  specific  internal  energy  of  propellant  at  the  exit  of  the  rocket  after  isentropic  expansion, 

absorption  volume. 

Vjfc*  normalized  absorption  volume,  = Vabs/(3. 

C*  normalized  coupling  coefficient,  C*  =  C/p. 
vi*  normalized  exit  velocity,  ve*  = 

(Uc-u0)*  normalized  specific  internal  energy,  (ue-u°)*  =  (u.-u0)/  J524>. 

T  temperature,  K. 

P  pressure,  bar. 

h  specific  enthalpy,  J/kg. 

s  specific  entropy,  J/kg  K. 

Mm  average  molecular  weight  of  a  mixture,  g/mole. 

Cp  specific  heat  capacity  at  constant  pressure. 
va  velocity  of  sound,  m/s. 

X(e-)  mole  fraction  of  electrons, 

s  expansion  ratio 

Aj  area  of  exit  surface 

At  area  of  sonic  surface  or  throat  area 

Subscripts,  Acronyms,  symbols 

i  initial  value  of  property. 

f  final  value  of  property. 

c  property  in  chamber 

t  property  in  throat 

e  property  in  exit  plane 

p  property  of  propellant 

<x>  denotes  mass  weighted  average  of  x 

HELSTF  High  Energy  Laser  System  Test  Facility 

PLVTS  Pulsed  Laser  Vulnerability  Test  System 

MLL  Myrabo  Laser  Lightcraft 

STP  Standard  Temperature  and  Pressure,  298  K,  1 .01326  bar 


INTRODUCTION 

Laser  propulsion  is  limited  by  laser  power,  so  optimization  of  the  laser  propulsion  mission  may  be  factored 
into  optimization  of  four  energy  conversion  efficiencies,  which,  in  a  first  approximation,  are  independent  of  each 
other.  In  this  idealization  the  kinetic  energy  of  the  propelled  vehicle  at  the  end  of  the  mission  may  be  expressed 
simply: 

(1)  Ef  =  '/zmfVf2  =  Ti<xpyEL. 

The  “propulsion  efficiency”,  t\,  is  the  efficiency  with  which  jet  kinetic  energy  is  converted  into  vehicle 
kinetic  energy.  Sutton1  pointed  out,  more  than  50  years  ago,  that  the  instantaneous  propulsion  efficiency  varies 
during  a  rocket  mission  and  that  it  is  unity  only  when  the  vehicle  velocity  in  the  inertial  frame  is  equal  to  the  jet 
velocity  m  the  rocket  frame.  Thus,  unit  propulsion  efficiency  is  achieved  only  when  the  jet  is  deposited  as  a 
stationary  mass  relative  to  an  observer  in  the  inertial  frame  of  reference. 

Then,  25-years  ago,  Moeckel2  and  Lo3  independently  and  nearly  simultaneously  published  analyses  of  the 
optimization  of  laser  rocket  propulsion  by  maximizing  the  overall  mission  average  efficiency  of  conversion  of  jet 


kinetic  to  vehicle  kinetic  energy.  Most  recently,  Phipps,  Reilly  and  Campbell  (2000, 200 1)4  cited  Moeckel’s  paper 
in  their  comprehensive  analysis  of  the  single  stage,  constant  I,p  Earth  to  LEO  rocket  mission.  They  reiterated  the 
fundamental  limit  that  Newton’s  second  law  imposes:  for  rocket  missions  that  start  at  zero  initial  velocity,  the 
maximum  q  is  0.648,  which  is  achieved  when  f  =  0.203  and  vf'v.  =  1.595.  For  the  Earth  to  LEO  mission  the 
effective  “delta  v”  (vf)  is  about  10  km/s,  so  the  optimum  single  stage  to  orbit  jet  velocity  is  -  6.27  km/s,  or  specific 
impulse  ~  640  s. 

In  this  paper  we  report  a  continuation  of  our  previous  work5  and  report  detailed  thermodynamics  analysis 
of  energy  conversion  in  the  “blowdown”  of  laser  heated  air  via  the  path  that  minimizes  entropy  increase.  The  results 
derive  from  analysis  of  isentropic  expansions  from  chemically  equilibrated  states  that  may  be  specified  by  their 
internal  energy,  density,  and  overall  stoichiometry.  Under  conditions  of  chemical  equilibrium  the  mixture 
composition  adjusts  itself  to  minimize  the  Gibbs  free  energy,  or  chemical  potential.  The  accuracy  of  the  equlibrium 
composition  depends  on  the  accuracy'  of  specification  of  relevant  species  in  the  mixture  and  their  thermodynamic 
properties.  The  analysis  spans  a  temperature  range  of  2000  to  24000  K,  and  several  orders  of  magnitude  in  pressure. 

We  discuss  and  analyze  measurements  of  the  overall  efficiency  of  conversion  of  laser  energy  to  propellant 
kinetic  energy,  a(3,  based  on  various  ballistic  pendulum  and  flight  experiments  with  Myrabo  Laser  Lightcraft,  MLL 
[Messitt,  Myrabo,  and  Mead  (2000)6;  Mead,  Squires,  Beairsto.  and  Thurston  (2000)'].  The  Phipps,  et  al.4  study 
defined  an  “ablation  efficiency  “  and  analyzed  the  Earth  to  LEO  mission  with  unit  ablation  efficiency.  Their 
ablation  efficiency  is  equivalent  to  the  product  of  our  a  and  p.  By  analyzing  experimental  results,  we  are  able  to 
narrow  the  range  of  p  that  operates  during  the  heating  process.  These  3  values  are  somewhat  larger  than  those 
reported  by  Wang,  et  al.8  for  CFD  plasma  models  of  the  heating  process.  It  has  been  pointed  out  that  P  approaches 
zero  as  the  plasma  temperature  approaches  ~  40,000  K,  where  the  plasma  frequency  approaches  the  laser 
frequency.8-  9 

The  Rocket  Equation 

The  thrust  that  results  from  expulsion  of  matter  at  velocity  ve  from  a  vehicle  of  mass  m  is  expressed  by 
Newton’s  second  law  as 


(2) 


F= 

dt 


where  mv,  is  the  momentum  of  the  jet  exhaust  in  the  vehicle  frame  of  reference,  [Corliss,  (I960)]’0.  For  the  case 
where  v*  is  constant, 


(3) 


F  = 


dm 


Equation  (2)  may  also  be  used  to  define  an  average  exit  velocity  for  rockets  where  ve  is  not  constant,  such 
as  blowdown  of  a  specified  mass  of  hot  propellant  from  a  fixed  volume,  e.g..  as  in  laser  rockets  and  pulse  detonation 
rockets: 


(4) 


<ve>  =  - 


t  mf  Pf 

jFdt  jd(mve)  jd(pve) 

0 _  mi  B _ 

mf  mf  Pf 

Jdm  jdm  Jdp 

m  j  nij  p. 


so  that  <ve>  is  the  mass  weighted  average  exit  velocity-.  Chemical  thermodynamics  may  be  used  to  rigorously 
establish  upper  limits  of  <ve>  when  the  propellant  equation  of  state  is  known  and  the  initial  and  final  states  of  the 
propellant  expansion  are  specified. 

The  Rocket  Equation  results  from  a  balance  of  the  force  exerted  by  the  propellant  on  the  vehicle  and  the 
motion  of  the  vehicle  under  the  influence  of  the  propulsive  force  as  required  by  Newton’s  second  law.  Thus,  in  the 
absence  of  other  forces,  such  as  body  (gravitational)  force  and  drag  force, 


4 


,,v  „  dm  dv 

<5)  F-V'irm7' 


where  v  is  the  vehicle  velocity  in  the  inertial  frame  of  reference,  i.e.,  the  velocity  relative  to  a  fixed  point  in  space. 
Elimination  of  time  in  Equation  (5)  yields  the  expression  for  conservation  of  momentum,  mdv  =  -vedm,  which  may 
be  integrated  between  the  limits  of  initial  and  final  mission  velocity  (v;  and  vf)  and  mass  (m*  and  mf)  to  produce  the 
Rocket  Equation, 

*  mf  T  vf "  vl 

(6)  f  =  —  =  exp  -  — — 1 

mj  L  ve 


Overall  Efficiency  of  Conversion  of  Laser  Energy  to  Propellant  Kinetic  Energy.  aB 

The  efficiency  of  conversion  of  laser  energy  to  propellant  kinetic  energy  may  be  defined  by  energy 
conservation  for  the  general  case  of  variable  ve  such  as  occurs  with  blowdown  expansion  of  laser  heated  air. 

(7)  Ep  =  =  o$  El, 

where  the  mass  weighted  average  of  the  square  of  the  propellant  exit  velocity  is 


Jfd(pve2) 


Pc 


The  impulse,  I  =  jFdt,  imparted  to  a  test  article  by  expansion  of  its  propellant  may  be  accurately  measured  with  a 
ballistic  pendulum.  Momentum  conservation  requires  equivalence  between  the  measured  impulse  and  the  propellant 
impulse  so  that 

(9)  I  =  nipCv^. 

The  momentum  coupling  coefficient,  also  a  measured  quantity,  is  the  impulse  imparted  to  a  test  article  per  unit  laser 
energy  incident  on  the  propellant, 

00)  0-1. 

el 

Using  the  definitions  embodied  in  Equations  (7)  -  (1 0),  C  may  be  expressed  in  terms  of  cc,  (3,  <v,>,  and  <ve2>: 

(id 

<ve>  L  <ve2  >  J  <ve> 

If  ve  is  constant,  <t>  =  <ve>2/<ve2>  =  1 .  Thermodynamics  may  be  used  to  rigorously  establish  inviolate  upper  limits 
to  <h  and  a  for  any  specified  free-expansion  blowdown  process  when  the  propellant  equation  of  state  is  known.  The 
<i>  factor  depends  on  the  mass  distribution  of  exit  velocities,  and  is  mathematically  limited  to  0.5  <  O  <  1.  It  will  be 
shown  (Figure  6,  vide  infra)  that  <J>  for  optimum  blowdown  of  laser  healed  air  to  1  bar  pressure  increases  from  0.95 
at  low  energy  (2  MJ/kg)  to  0.98  at  high  energy  (60  MJ/kg).  The  $  factor  arises  in  Equation  (11)  because  the 
measured  quantity,  the  jet  impulse,  is  proportional  to  mass  weighted  average  velocity  whereas  the  jet  kinetic  energy 
is  proportional  to  the  mass  weighted  average  of  the  squared  velocity. 
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Experimental  determination  of  aB 


The  value  of  a(3  may  be  determined  within  a  factor  of  <J>  <  1  by  experimental  measurement  of  the  impulse 
imparted  to  a  test  article  (I)  when  a  laser  pulse  of  known  energy  (EL)  ejects  a  known  amount  of  propellant  mass 
K): 


(12) 


apd)  = - 

'  2mpEL 


Cl  C<ve> 
2mp“  2 


l<ve> 

2El 


Thus,  measurement  of  cep  requires  knowledge  of  the  propellant  mass  that  is  associated  with  the  measured 
impulse  and  laser  pulse  energy.  In  the  absence  of  a  mass  measurement,  a  low'er  limit  to  propellant  mass  and  an 
upper  limit  to  exit  velocity  may  be  established.  Since  P  <  1  is  required  for  energy  conservation, 


(13) 

r 

m  >  — * — 
p  2a4>EL 

(14) 

2a*EL 
ve  S  ,  • 

and 


Use  of  thermodynamics  to  establish  upper  limits  to  <t>  and  a  enables  additional  restrictions  to  be  placed  on  the 
permissible  experimental  values  of  the  upper  limit  vt  and  lower  limit  m?. 

Maximum.  Efficiency  of  Conversion  of  Propellant  Internal  Energy'  to  Propellant  Kinetic  Energy,  a. 

In  the  analysis  of  idealized  thermodynamic  expansion  of  laser  heated  air,  the  notion  of  an  energy  absorption 
volume  is  invoked  that  contains  a  mass  of  propellant  mp  =  pcV*,  into  w’hich  an  amount  of  energy  PEL  is  deposited. 
The  time  scale  for  energy  absorption  is  much  shorter  than  that  for  expansion  so  that  the  propellant  density  within 
VabS  (the  chamber)  remains  constant  during  energy  absorption.  This  enables  the  initial  specific  internal  energy  of  the 
propellant  to  be  defined, 

(15)  uc  -  u°  =  PEL/pcVlbs, 

where  pc=  1.18kg/m3  andu°  =  0.09x  106  J/kg  for  air  at  STP.  Table  I  providesa  convenient  list  of  values  of  the 
normalized  absorption  volume,  V^*  =  Vlbs/p,  derived  from  Equation  (15)  for  various  values  of  14  -  u°and  EL. 

Figure  1  shows  a  cross-section  of  the  test  article  with  a  ring  of  Delrin  installed  in  the  shroud.  The  Delrin 
shown  occupies  a  volume  of  7  cmJ,  which  may  be  used  to  visualize  a  reasonable  absorption  volume  for  the  case 
where  Delrin  is  absent  and  air  is  the  heated  material.  Table  1  shows  that  a  similar  absorption  volume  for  air  would 
produce,  with  unit  p  and  nominal  El  values  between  100  and  400  J,  heated  air  with  1 0  to  40  MJ/kg  internal  energy. 

If  the  Delrin  surface  shown  in  the  figure,  about  25  cm3,  is  a  suitable  representation  of  the  sonic  surface  of  expanding 
air,  then,  with  an  idealized  plug-nozzle  exit  area11  of  -  350  cm3,  the  expansion  ratio  in  this  test  article  may  be  as 
large  as  ~  14. 

Perfect  isentropic  conversion  of  internal  energy’  to  propellant  kinetic  energy  occurs  with  no  losses  so  that 

(16)  <ve2>=  2<uc  -  ut>  =  2a(uc  -  u°),  where 

(17)  a  =  <uc-ue>/(uc  -  u°)  =  <v,3>/2(uc  -  u°),  and 

(18)  <ve>  =  2<(uc-uc),->. 

These  definitions  generate  a  second  expression  for  C  in  terms  of  the  specific  energy'  of  laser  heated  propellant: 

(19)  C  =  P[2a<J>/(uc-u°)],/3. 
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Figure  2  shows  the  relationships  between  six  variables  of  interest:  C,  a,  p,  $,  <ve>,  and  [uc  -  u  ]  that  derive  from 
definitions  given  by  Equations  (1 1)  and  (19).  The  Figure  shows  C*  =  Cl  p  vs  a  plot  with  lines  of  constant  <ve>*  = 
<ve>/®  and  lines  of  constant  (uc  -  u0)*  =  (ue  -  u°)/<D.  As  one  proceeds  from  die  origin  along  a  paraboloidal  line  of 
constant  (uc  -  u0)*,  which  is  also  a  line  of  constant  entropy,  both  C*  and  a  increase.  At  constant  a,  C*  decreases  as 
<ve>*  and  -  u0)*  increase.  Knowledge  of  C  and  a  fixes  values  of  <Vj>/  (3$  and  [ue  ■  u°]/p"d>.  Figure  1  may  also 
be  interpreted  as  a  C  vs  a  plot  with  lines  of  constant  <ve>/  |3<t>  and  lines  of  constant  [uc  -  u°]/f32<E>.  The  absorption 
volume  may  also  be  factored  into  this  parameter  space  with  use  of  Equation  (15)  or  Table  1 . 

Thermodynamic  Limitations  to  a  and  d>. 

Figure  3  shows  the  chemical  equilibrium  Mollier  diagram  for  air  up  to  24,000  K.  Fipre  3  is  based  on  the 
database  maintained  at  NASA/Glenn  [McBride  and  Gordon  (1996)]12,  which  is  certified  accurate  up  to  20,000  K 
and  which  is  based  on  extended  9-parameter  fits  to  enthalpy,  heat  capacity,  and  entropy  of  neutral  species  and  singly 
charged  ions.  Above  20,000  K  doubly  charged  ions  begin  to  contribute  but  these  are  not  included  in  the  database. 
This  limitation  leads  to  predictions  of  temperatures  (at  specified  u  and  p)  that  are  too  high  for  plasmas  above  - 
20,000  K. 

Figure  4  shows  a  series  of  vertical  lines  on  the  Mollier  diagram.  These  are  representations  of  equilibrium 
isentropic  expansions  that  originate  from  initial  states  located  along  the  constant  density  line,  p  =  1.18  kg/m5,  and 
specific  internal  energies  ranging  from  1  to  100  MJ/kg.  Table  2  summarizes  other  thermodynamic  properties  of 
these  initial  equilibrium  states  of  interest:  T,  P,  h,  s,  Mm,  Cp,  v„  Cp/Cv,  and  X_(e").  Table  3  provides  a  similar  summary 
of  properties  for  the  case  of  Mach  5  air  at  its  stagnation  density1,5,  5.9  kg/m5.  Since  the  entropy  of  the  initial  and 
final  states  are  equal,  the  thermodynamic  state  of  the  propellant  in  the  exit  surface  is  uniquely  defined  when  only 
one  additional  property  in  the  exit  surface  is  specified,  such  as  the  exit  pressure  or  the  expansion  ratio,  which  are 
indicated  in  Figure  4.  The  expansion  ratio,  e,  is  the  ratio  of  the  area  of  the  exit  surface  to  the  area  of  the  sonic 
surface  or  nozzle  throat,  and  for  isentropic  expansions  this  may  be  represented  in  terms  of  thermodynamic  properties 
in  the  nozzle  throat  and  exit  plane:  e  =  Ae/A,=  p,Vpeve . 

Optimum  Blowdown  to  P.  =  1  bar.  Mass  Weighted  Average  Quantities 

Figure  5  shows  a  representation  of  blowdown  from  the  initial  state  where  u^u,,  =  2E3  kJ/kg  and  die  initial 
density  is  that  of  STP  air.  The  series  of  vertical  lines  are  located  at  equally  spaced  density  increments.  The 
instantaneous  and  mass  weighted  average  quantities  based  on  Equations  (1 6)  -  (19)  are  shown  in  Figures  6-9. 

Figures  10  and  1 1  show  the  transformations  of  the  isentropes  in  the  Mollier  plane  (u-s  plane)  to  the  C*-a 
plane.  Lines  of  constant  s  and  pe  are  almost  exactly  coincident.  Lines  of  constant  exit  pressure  run  nearly  parallel  to 
lines  of  constant  e  and  p*  and  all  are  nearly  vertical,  indicating  that  alpha  is  nearly  independent  of  ve  and  u^-u0. 

DISCUSSION 


Coupling  coefficients  measured  with  the  Figure  1  test  model  were  reported  in  our  previous  paper5.  With 
increasing  laser  energy  they  rise  to  plateaus  above  about  300  J.  At  El  ~  300  J,  C(Delrin)  ~  350  Ns/MJ  and  die 
ablated/vaporized  mass  was  mp  ~  35  mg.  This  means  that  <v,>  ~  3000  m/s  by  Equations  (9)  and  (10),  and  aP<J>  ~ 
0.5  by  Equation  (1 1).  Thus,  a  >  0.5,  which  is  remarkably  high.  Air  and  Delrin  will  show  veiy  similar  expansion 
behavior.  As  shown  by  Figure  1 1,  the  dependence  of  a  on  the  density  of  the  heated  air  is  quite  strong.  At  [Uc-u0]*  = 
10  MJ/kg  and  expansion  to  1  bar,  Figure  1 1  shows  that  the  instantaneous  a  increases  from  about  0.43_to  about  0.60 
when  the  density  increases  from  the  STP  value  (1.1 8  kg/m3)  to  its  Mach  5  stagnation  value  (5.9  kg/m3).  The 
instantaneous  a  values  decrease  to  about  032  and  0.5  respectively  when  the  mass  weighted  average  a  for  the  free- 
expansion  blowdown  process  is  calculated^  An  a  value  around  03  is  reasonable  when  die  density  of  the  ablated 
and  vaporized  delrin  is  as  high  as  ~  6  kg/m3  and  die  blowdown  expansion  is  near  perfect  Most  importantly,  it 
appears  that  most  of  the  inefficiency  in  the  composite  efficiency  is  carried  by  a  and  that  f3<t>  is  very  close  to 
unity. 

Coupling  coefficients  at  EL  >  300  J  for  air  were  found  to  depend  strongly  on  the  quality  of  the  laser  beam, 
as  between  a  tightly  focused  beam  that  produced  lower  C  ~  1 00  Ns/MJ  than  a  loosely  focused  beam,  which 
produced  C  ~  150  Ns/MJ.  This  may  be  due  to  the  tight  beam  heating  a  smaller  mass  of  air  to  a  higher  energy  than 
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Ae  more  diffuse  loosely  focused  beam.  Although  the  exit  velocity  would  be  higher  in  the  tight  beam  case,  the  total 
impulse  may  be  lower  because  the  heated  mass  is  lower.  Figure  1  shows  the  geometry  and  size  relationship  of  a  7 
cm3  absorption  volume  inside  the  shroud,  which  contains  ~  8  mg  of  air.  With  C(air,  loose  focus)  =  150  Ns/MJ  at  EL 
=  300  J,  and  mp  =  8  mg  we  may  deduce  <vt>  ~  5600  m/s,  and  a30  ~  0.42.  If  the  absorption  volume  is  double,  then 
<ve>  and  apd>  are  halved.  Figure  10  shows  that  air  heated  to  10  MJ/kg  for  example  would  blowdown  to  1  bar  with 
a  =  0.32  (equilibrium  expansion)  or  a  =  0.27  (frozen  expansion). 

If  we  accept  that  a  reasonable  upper  limit  operational  alpha  is  -  0.30  in  our  experiments  a  <  0.3,  then  the 
measured  C~  150  Ns/MJ  and  Equation  (9-11)  with  (30  =  1  require  <ve>  <  4000  m/s,  and  mp  >  11  mg.  Now  if  pq>  is 
~  0.3  as  has  been  suggested  by  CFD  modeling,8  then  the  upper  limit  of  <vc>  decreases  to  1200  m/s  and  the  lower 
limit  of  mp  increases  to  36  mg.  It  would  seem  apparent  that  the  value  of  (5  is  somewhat  larger  than  0.3  because  both 
the  upper  limit  <ve>  and  lower  limit  mp  are  not  reasonable  for  the  geometry-  shown  in  Figure  1 . 

CONCLUSIONS 


Experimental  studies  of  the  200-3/4  mode!  Myrabo  Laser  Lightcraft  with  air  heated  by  1 0.6  u  radiation 
from  a  C02  laser  showed  that  energy  conversion  efficiencies  of  laser  energy  to  propellant  kinetic  energy  were  at 
least  30%.  This  was  found  to  be  consistent  with  highly  simplified  analy  sis  of  equilibrium  (isentropic)  expansions 
from  initial  states  that  are  specified  by  a  single  parameter,  the  volume  into  which  the  laser  energy  is  absorbed.  The 
measured  exit  velocity  based  on  estimated  air  mass  is  in  the  neighborhood  of  3000  m/s.  It  should  be  noted  that 
beam  quality  plays  an  important  role  in  the  performance  of  the  model  200-3/4  MLL  which  is  counterintuitive 
inasmuch  as  lower  beam  quality  (energy  spread  over  a  larger  area)  produces  higher  C  coupling  coefficients. 
Computational  fluid  dynamics  modeling  of  the  absorption  (and  reflection)  of  laser  energy  and  expansion  of  the 
formed  plasma8  have  recently  been  carried  out.  The  simple  analysis  presented  here  may  only  be  useful  in  providing 
upper  limitations  to  the  conversion  of  laser  energy'  to  propellant  kinetic  energy  and  to  provide  a  simplified 
visualization  and  description  of  the  processes  occurring  in  blowdown  of  laser  heated  propellants  in  MLL  devices. 

Expansion  of  a  propellant  mass  that  was  heated  at  constant  volume  was  examined  under  conditions  where 
either  chemical  equilibrium  or  frozen  composition  tvas  maintained.  For  expansion  with  an  effective  area  ratio  of  ~ 
4,  which  is  appropriate  for  the  MLL,  a  maximum  of  25  to  50%  of  the  internal  energy  is  predicted  to  be  convertible 
to  propellant  kinetic  energy,  based  on  the  minimization  of  the  entropy  gain  of  the  blowdown  process.  With  the 
small  effective  area  ratio  -  4,  equilibrium  expansion  was  only  slightly  more  efficient  than  frozen  expansion. 

Heating  of  propellant  to  highly  ionized  states  resulted  in  lower  efficiency  energy  conversion  but  higher  exit  velocity. 
The  thermodynamic  limitations  are  illustrated  by  process  representations  of  blowdown  in  the  Mollier  plane. 
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Table  1.  Normalized  absorption  volume  for  air  at  1.18  kg/nv 
function  of  internal  energy  and  laser  energy. 
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a  surface  area  ~  25  The  iilcatizcil plug  nozzle  exit  area  is  ~  350  cm3. 


Table  2.  Thermodynamic  properties  of  equilibrium  air,  p  =  1.18  kg/mJ. 
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Table  3. 

Thermodynamic  properties  of  Mach  5  air  at  stagnation  density,  p 

=  5.90  kg/m3 
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Figure  2.  Defined  relationships  between  six  variables  of  interest:  C*.  a,  f),  (I>,< ve>  *.  and  /nc  -  u"J *  The  plots  show  C*  =  C/f  as  a 
function  of  a,  with  Hues  of  constant  vf  =  <veW<Pand  constant  /uc-u"f*  =  fu,-u"J/d>.  The  plots  may  alternatively  he  interpreted  as  a  C 
u*  a  plots  with  lines  of  constant  ve*  =  <v ,>/(]<!>  and  constant  luc-n"l*  =  hc-uV/  (ftp. 


specific  entropy  J/kg  K 

Figure  3.  MoUier  diagram  for  air  including  singly  ionized  species.  Molecular  weights  are  inMcated  at  intersections  of 
isobars  and  isotherms.  The  lower  diagram  shows  a  heavy  constant  density  line,  p  - 1.18  kg/m3  above  a  heavy  constant 
pressure  line,  P-1  atm.  The  maximum  energy  initial  states  of  laser  heated  STP  air  lie  on  the  constant  density  line  and 
the  optimally  expanded  states  lie  vertically  below  on  the  constant  pressure  line 


phi  =  <ve>z/<vez>.  exitvebcity  (mfc)  alpha 


Instantaneous  Quantities 


Figure  6.  Instantaneous  quantities  for  equilibrium  blowdown  of  heated  air from  initial 
density  of  1.18  kg/m3  and  specific  internal  energies  ranging  from  2  to  60  MJ/kg  to  a  final 
pressure  of  1  bar. 


cou  pling  CO  efficient  ( N  s/MJ )  exit  v  eto  city  ( mfe ) 
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density  (kg/in3) 

Figure  7.  Mass  weighted  average  quantities  for  equilibrium  blowdown  of  heated  air  from 
initial  density  of  1.18  kg/m3  and  specific  internal  energies  ranging  from  2  to  60  MJ/kg  to 
a  final  pressure  of  1  bar. 
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Figure  9.  Comparison  of  mass  weighted  average  FROZEN  and  EQUILIBRIUM 
blowdown  quantities  from  2  and  10  MJ/kg  initial  internal  energy.  Initial  density  US 
kg/m3  with  blowdown  to  J  bar  external  pressure 


obtained  from  initial  {uc~h°/*  states  of  2E3,  and  I  Ed  J/kgfor  the  frozen  expansion  and  2E3,  6E3, 1E4,  and  4E4  kJ/kg  for  the  equilibrium  expansion. 
The  results  of  the  two  frozen  blowdown  Integrations  to  Pa„  - 1  bar  are  plotted  with  those  of  the  equilibrium  blowdown  to  show  that  the  differences  in 
alpha  are  small,  l.e.,  at  low  energy  (2E3)  0.30  and  0.29  and  at  high  energy  (1E4)  0.32  and  0.27 for  equilibrium  and  frozen  blowdown,  respectively. 


rO 

tO  <s 

ui 

< 

< 

ui 

cP 

.  .c?  f 

^  <s 

to  < 


*1^  ^  ""  c—  ^  *5  ■  .  ^J:.  '?  a 

\.  \ :  \” “ “ r « »  ’• .  '*  «vvv{^5|).-A 

.  r*\*'-V.'v  •  ■  'vs.  ^  •'••>  V  -f 

Ai  ~  °  Aj'  ’  *P  ;.  ’*%  .  ^  .  ;  '7^^  '  ■- 

H -:V  v ; ; -'-' ;•$'  v"' \f^v*&9a * ’ 

A '  :V- v  V '  I  ^  ~  ^ '  'V  ‘ ■'  ’  S3 1 

.  • .' \  *  v°:93 : 

Cv'r- 


V)  'r , 

ra  ^r>. 

a  c/. 

x  . 
a 

E 

•2  „  .'■ 

5 


^TX~^r,?^rvrrr 

'V  •  \  •  \  '  “s'~'o.\  '.  ' 

■'. ;  \  \ : -\  v*rv 


—  - ... .  .N„ ' 


%  *  ❖ 


rifl/SN  'e}9q/o=.o  ‘juaa^soo  6in|dnoo  P92|elujou 

•C  \p 

1  ^  til  u:  t&uiu)  STC>i/r>)).(cr|-Sr') 

r-  <-  c4  <■  iD  tD  r-  CP 


(P  a>  vi.  •  ,  a 

.T-.  Vw  5>  ••'  ■  :  f<  {S/Ui»).  A 


^  \  *,'w  •'  •••'  '.  ••  <  (S/UJ))) .  A 

N.-.  •>-»'■:  •. :- ^ 

r-^=~ .  \  ».A  ”.  i'i.  "'a.  -  .:. 

V  '<••  ?'<-  ^*<83i 

-  „  .  . ... 

®  °c  ‘V.  V V .  Y~-» !' '  '' •' '■-  ■  v:  v  ■  ->zo i 

^  ^  ■  '^^rv^Cv4.-''''v>- 

''  •  >•,.  '-c  ."4  q's"  ••i«9S2 

\;  \  ‘  '••'  *'%-« -  j  \.  ' '*  C>'  Vv  V-  e*'~ZZ  O 
Jeq  Z  =  ’d 


S  ^ 

a 

x 

<D 

E 

•B 

n 


riAl/SN  'e}9q/0=O  'jusoiiEoo  6ui|dnoo  psaieiujou 


Figure  1 1.  Comparison  of  Equilibrium  expansion  from  laser  healed  STP  air  (I.  IS  kg/m')  and  Mach  5  air  al  stagnation  density  (5.9  kg/m').  In  the 
SIT  air  diagram  (on  left),  the  circles  and  nearby  crosses  represent  the  blowdown  quantities  obtained  from  initial  fite-u"J*  states  of2EJ,  and  !E4  J/kg 
for  the  frozen  expansion  and  2E3.6E3, 1E4,  and  4 E4  kJ/kg  for  the  equilibrium  expansion. 


